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The interaction of polymers with turbulent shear flows is examined. We focus on the structure 
of the elastic stress tensor, which is proportional to the polymer conformation tensor. We examine 
this object in turbulent flows of increasing complexity. First is isotropic turbulence, then anisotropic 
(but homogenous) shear turbulence and flnally wall bounded turbulence. The main result of this 
paper is that for all these flows the polymer stress tensor attains a universal structure in the limit 
of large Deborah number Ife S> 1. We present analytic results for the suppression of the coil-stretch 
transition at large Deborah numbers. Above the transition the turbulent velocity fluctuations are 
strongly correlated with the polymer's elongation: there appear high-quality "hydro-elastic" waves 
in which turbulent kinetic energy turns into polymer potential energy and vice versa. These waves 
determine the trace of the elastic stress tensor but practically do not modify its universal structure. 
We demonstrate that the influence of the polymers on the balance of energy and momentum can 
be accurately described by an effective polymer viscosity that is proportional to to the cross-stream 
component of the elastic stress tensor. This component is smaller than the stream- wise component 
by a factor proportional to Ite^. Finally we tie our results to wall bounded turbulence and clarify 
some puzzling facts observed in the problem of drag reduction by polymers. 
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I. INTRODUCTION 

The dynamics of dilute polymers in turbulent flows is 
a rich subject, combining the complexities of polymer 
physics and of turbulence. Besides fundamental ques- 
tions there is a significant practical interest in the sub- 
ject, particularly because of the dramatic effect of drag 
reduction in wall bounded turbulent flows [1-4]. On the 
one hand the polymer additives provide a channel of dis- 
sipation in addition to the Newtonian viscosity; this had 
been stressed in the past mainly by Lumley [5, 6]. On 
the other hand the polymers can store energy in the form 
of elastic energy; this aspect had been stressed for exam- 
ple by Tabor and De Gennes [7] . A better understanding 
of the relative roles of these aspects requires a detailed 
analysis of the dynamics of the "complex fluid" obtained 
with dilute polymers in the turbulent flow regime. 

Some important progress in the theoretical description 
of the statistics of polymer stretching in homogeneous, 
isotropic turbulence of dilute polymer solutions was of- 
fered in Refs. [8-10]. Here we want to stress the fact 
that in the practically interesting turbulent regimes, and 
in particular when there exists a large drag reduction ef- 
fect, the characteristic mean velocity gradient (say, the 
mean shear, Sq), is much larger than the inverse polymer 
relaxation time, Tp, SoTp ^ 1. Usually this parameter is 
referred to as the Deborah or Weissenberg number: 



(1.1) 



Indeed, the onset of drag reduction corresponds to 
Reynolds number Tiecr at which De 1. Large drag 
reduction corresponds to 7?e 3> Tiecr at which Tk ~ 
'Re/'Re„ > 1. 

The aim of this paper is to provide a theory of the poly- 
mer stress tensor 11 [defined in Eq. (2.1b)] in turbulent 
fiows in which Tb ^ 1. The main result is a relatively 
simple form (1.2a) of the mean polymer stress tensor Ho 
that enjoys a high degree of universality. Denoting by 



X and y the unit vectors in the mean velocity and mean 
velocity gradient directions (streamwise and cross-stream 
directions in a channel geometry) and by z = x x y the 
span-wise direction, we show that tensor Ho has a uni- 
versal form: 
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fori:te>l. (1.2a) 



One sees that in the (x-y)-plane the tensorial structure 
of IIo is independent of the statistics of turbulence: 



(1.2b) 



Due to the to the symmetry of the considered flows with 
respect to reflection z ^ ~z, the off-diagonal compo- 
vanish identically: 



nents Ilg^ and IIq^ 



JTXZ 
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(1.2c) 



The only non-universal entry in Eq. (1.2a) is the dimcn- 
sionless coefficient C. We show that C = 1 for shear 
fiows in which the extension of the polymers is caused by 
temperature fluctuations and/or by isotropic turbulence. 
For anisotropic turbulence the constant C is of the order 
of unity. 

The universal tensorial structure (1.2) has important 
consequences for the problem of drag reduction in dilute 
polymeric solutions. We show that the effect of the poly- 
mers on the balance of energy and mechanical momentum 
can be described by an effective polymer viscosity that is 
proportional to the cross-stream component of the elastic 
stress tensor Hq^. According to Eq. (1.2b) this compo- 
nent is smaller than the stream-wise component Ilg^ by 
a factor of 2Xte^. This flnding resolves the puzzle of the 
linear increase of the effective polymeric viscosity with 
the distance from the wall, while the total elongation of 
the polymers [dominated by the largest component Ilg^ 
of the tensor (1.2a)] decreases with the distance from the 
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wall. Wc show that our results for the profiles of the 
components IIq-' are in a good agreement with DNS of 
the FENE-P model in turbulent channel flows. 

The paper is organized as follows. In Sec. II A we 
present the basic equations of motion of the problem. 
In Sec. II B we describe the standard Reynolds decom- 
position of the equations of motion for the mean and 
fluctuations of the relevant variables. In Sec. II C we find 
general solution for IIo and formulate (the very weak) 
conditions at which it simplifies to the universal form 
(1.2). 

In Section III we analyze the phenomenon of polymer 
stretching in shear flows with prescribed turbulent ve- 
locity fields, cither isotropic, Sec. IIIB, or anisotropic. 
Sec. inc. We relate the value of C to the level of 
anisotropy in the turbulent flow. We also show in this 
section that a strong mean shear suppresses the thresh- 
old value for the coil-stretched transition by a factor Ve^, 
see Eqs. (3.10). 

Sections IV and V are devoted to the case of strong tur- 
bulent fluctuations for which the fluctuating velocity and 
polymeric fields are strongly correlated. In Sec. IV we 
consider the case of homogeneous shear flows. In Sec. V 
we discuss wall bounded turbulence in a channel geom- 
etry, compare our results with available DNS data and 
apply our finding to the problem of drag reduction by 
polymeric solutions. 

Section VI presents summary and a discussion of the 
obtained results . 

Appendix A offers some useful exact relationships for 
the energy balance in the system. 

II. EQUATIONS OF MOTION AND SOLUTION 
FOR SIMPLE FLOWS 

A. Equations of motion for dilute polymer 
solutions 



a single relaxation time Tp, which is constant at small 
extensions. In this case 11 can be written as 



n(t, r) = -^R{t, r)R{t, r) , (2.1b) 

where z^op is polymeric viscosity for infinitesimal shear, 
and R is the polymer end-to-end elongation vector, nor- 
malized by its equilibrium value (in equilibrium RR = 
S). The average in Eq. (2.1b) is over the Gaussian statis- 
tics of the Langevin random force which describes the 
interaction of the polymer molecules with the solvent 
molecules at a given temprature, but not over the tur- 
bulent ensemble. 

b. The equation for the elastic stress tensor U has 
the form 

^ + (V • v)n = s • n + n • - 1 (n - u,^) , 

dt Tp 

S = {Wy , S'^ ^ dV'/dx^ , (2.1c) 

Ileq = Ileq^ , Ileq = fQp/Tp . 

Here S = S{t, r) is the velocity gradient tensor and S is 
unit tensor. 

c. Choice of coordinates. In this paper we consider 
both homogeneous and wall bounded shear fiows. In both 
cases we choose the coordinates such that the mean ve- 
locity is in the x direction and the gradient is in the y 

direction. In wall bounded flows x and y are unit vectors 
in the streamwise and the wall normal directions: 

Vo = {V)=So-r, 5o = 5oxy, (2.2) 
only y^f = 5oy ^ and S^"" = So ^ . 

The unit vector orthogonal to x and y (the span-wise 
direction in wall bounded cases) is denoted as z. 



a. The equation for the velocity field of the complex 
fluid reads: 

^ + {V ■V)V = i/qAV- VP-h v-n, 

V V = . (2.1a) 

Here V = V{t,r), P = P{t,r) and vq are the fluid ve- 
locity, the pressure, and Newtonian viscosity of the neat 
fluid respectively. The fiuid is considered incompressible, 
and units are chosen such that the density is 1. The ef- 
fect of the added polymer appears in Eq. (2.1a) via the 
elastic stress tensor 11 = n(t, r). 

In this paper we consider dilute polymers in the limit 
that their extension due to interaction with the fluid 
is small compared with their chemical length. In this 
regime one can safely assume that the polymer extension 
is proportional to the applied force (Hook's law). We 
simplify the description of the polymer dynamics assum- 
ing that the relaxation to equilibrium is characterized by 



B. Reynolds decomposition of the basic equations 

In the following we need to consider separately the 
mean values and the fluctuating parts of the velocity, 
V{r,t), the velocity gradient S{r,t), and the elastic 
stress tensor Il{r,t) fields. We define the fluctuating 
parts via 

V = Vo+v, S = So + s, n = no + n,... (2.3) 

All the mean quantities will be denoted with a "o" sub- 
script (e.g. the mean velocity Vq), and all the fluctuating 
quantities by the lower-case letters v, p, tt, s, etc. These 
mean values are computed with respect to the appropri- 
ate turbulent ensemble. Note that the mean pressure 
in a homogeneous shear flow is zero, Pq = 0, and all 
mean quantities (except the mean velocity Vq = Soy, 
of course) are coordinate independent and, by definition, 
time independent. 
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1. Equations for the mean objects 

a. The equation for the mean elastic stress tensor 
Ho follows from Eq. (2.1c): 



^ + -]Uo = So-Uo + no-S^ + Q, (2.4a) 



Q = — Heq + Qt , 
Qt = (s • TT + TT • s"^) , 



(2.4b) 
(2.4c) 



where Heq is given by Eq. (2.1c), and (• • • ) denotes an 
average with respect to the turbulent ensemble. In com- 
ponents: 

+ ^ ) = + + ' 

= — Heq^'^' + (s''=7r*=^' + s^'^tt''') . (2.4d) 
Here D/Dt is the mean substantial time derivative 

The substantial derivative vanishes in the stationary case, 
when all the statistical objects are t- and .T-indcpcndcnt. 

b. The equation for the mean velocity follows from 
Eq. (2.1a): 

^ = ^[,,So-w-y + iiiy], (2.6) 

were W is the Reynolds stress tensor: 

W'^ = (v'v^) . (2.7) 

In the stationary case dVo/dt = and Eq. (2.6) gives 

i^oSo - + = P , (2.8) 

where constant of integration V has a physical meaning 
of the total momentum flux. The three terms on the LHS 
of Eq. (2.8) describe the viscous, inertial and polymeric 
contributions to V. In a homogeneous-shear geometry 
all theses terms are ^-independent constants, whereas in 
a channel geometry they depend on the distance to the 
wall. Notice, that Eq. (2.8) also can be considered as an 
equation for the balance of mechanical forces in the flow. 



2. Equations for the fluctuations 

The equations for the fluctuating parts v and tt read: 
Dv 



Dt 

D 1 



So ■ V + uqAv — Vp -|- V • (tt — vv) , 

(2.9a) 

TT = So ■ n + IT ■ + s ■ Uo + no ■ s'^ 

+ , (2.9b) 



where V • = 0, and 

= {s ■ TT + TT ■ S"^) - V ■ VTT . (2.9c) 

Here (...) denotes the "fluctuating part of (...). 

C. Solutions for simple flow conflgurations 

1. Implicit solution for shear flows 

Note that the mean shear tensor Sq satisfies, besides 
the incompressibility condition TrjS'o} = Sq = 0, one 
additional constraint 



Q. 



Jk _ 







(2.10) 



Remarkably, this property of Sq uniquely distinguishes 
shear flows from other possibilities (elongational or ro- 
tational flows): if (2.10) holds, one can always choose 
coordinates such that the only nonzero component of So 

Sq. 

Having Eq. (2.10) consider Eq. (2.4a) in the stationary 
state: 



is s^,y 



Tp [So Uo + Uo S^ + Q) 



(2.11) 



We can proceed to solve this equation implicitly, treat- 
ing Q on the RHS as a given tensor, and solve the lin- 
ear set of equation for Ho. In the considered geometry 
Eq. (2.11) is a system of 4 linear equations, and the so- 
lution is expected to be quite cumbersome. However, 
property (2.10) allows a very elegant solution of this sys- 
tem. Using (2.11) iteratively (i.e., substituting instead of 
Ho on the RHS of Eq. (2.11) the whole RHS), we get 



Ho = STpS'o • IIo • 



tI {So-Q + Q- S^) + TpQ 



(s'o-Q + Q-{s'oy) 



One sees that due to Eq. (2.10) the last term in this 
equation vanishes. Repeating this procedure once again, 
we obtain the exact solution of (2.11) in the form of a 
finite (quadratic) polynomial of the tensor So- 

Ho = 2 Tp3 So • Q • S J + Tp2 ( So • Q + Q • So" ) + Tp Q . 

(2.12a) 

The individual components of the solution Eq. (2.12a) 
are given by: 



'-'-0 

TT^y 

TTXZ 

^^0 



= Tp(2Tk'Qy 



TpiVeQ 



vv 



^vv 



+ 2VeQ''y 



(2.12b) 



yrzz 



yiV^ 



. 



Notice that the xz and yz components of Ho vanish due 
to the symmetry of reflection z ^ —z, which remains 
relevant in all the flow configurations addressed in this 
paper. 
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In the Hmit Xte ^ 1 the tensorial structure (2.12b) can 
be simphfied to the form (1.2) if 

Q'^y <ii.VeQyy , Q'''' -^Ve^Qyy . (2.13) 

Indeed, in this case one can neglect the last two terms 
in IIq^ and the second term in IIq^ and obtain Eq. (1.2) 

with c = Q'^'/Qyy. 



In the following Sects. Ill and IV we will find explicit 
solutions for the clastic stress tensor in the presence of 
turbulence which share a structure similar to Eq. (2.15). 
Namely, for Tt ^ 1 the leading contribution to each 
component of Ho can be presented as Eq. (1.2), in which 
C, given by Eqs. (3.13) and (5.8), is some constant of the 
order of unity, depending on the anisotropy of turbulent 
statistics. 



2. Explicit solution for laminar homogeneous shear flows 

The solution (2.12) allows further simplification in the 

case of laminar shear flow, in which Qt of Eq. (2.4c) 
vanishes, Qt = 0, and thus, according to Eqs. (2.4b) and 
(2.1c), the tensor Q becomes proportional to the unit 
tensor: 



With this relationship Eq. (2.12a) simplifies to 



Ho 




(2.14) 



(2.15) 



Equation (2.15), in which Ileq is given by Eq. (2.1c). 
is important in itself as an explicit solution for IIo in the 
case of laminar shear flow. But more importantly, this 
equation together with Eq. (2.12) gives a hint regard- 
ing the tensorial structure of IIq even in the presence 
of turbulence, when the coupling between the velocity 
and the polymeric elongation field, leading to the cross- 
correlation Qt, cannot be neglected. 

To see why the simple result (2.15) may be relevant 
also for the turbulent case, note first that the non- 
diagonal elements with z projection, i.e. IIq^ and IIq^, 
are identically zero in general as long as the symmetry 
z ^ —z prevails. Second, for large Deborah numbers, 
rte ^ 1, the nonzero components in both Eqs. (2.15) 
and (2.12b) have three different orders of magnitude: 
n§^ > Hq^ > ~ in other words, the tensor Hq 
is strongly anisotropic, reflecting a strong preferential ori- 
entation of the stretched polymers along the stream- wise 
direction x. The characteristic deviation angle (from the 
X direction) is of order 0(1 /Ife). 

Notice that Eqs. (2.12) relate the polymeric stress ten- 
sor Ho to the cross-correlation tensor Qt, Eq. (2.4c). In 
its turn this tensor depends on the polymeric stretch- 
ing, which is described by the same tensor Ho. There- 
fore, generally speaking, Eqs. (2.12) remains an im- 
plicit solution of the problem, which in general requires 
considerable further analysis. However, as we discussed 
above, Eq. (2.12a) is more transparent than the start- 
ing Eq. (2.4). In particular, if the tensor Q is not 
strongly anisotropic, the tensorial structure of IIo is close 
to Eq. (2.15) for the laminar case. We will see below that 
this structure is indeed recovered under more general con- 
ditions. 



III. POLYMER STRETCHING IN THE 
PASSIVE REGIME 

A. Cross-correlation tensor Qt in Gaussian 
turbulence 

When Xfe ^ 1 the characteristic decorrelation time of 
turbulent fluctuations, Tcor ^ l/5'o, is much smaller than 
the polymer relaxation time Tp. Accordingly the turbu- 
lent fluctuations can be taken as ^-correlated in time: 



(3.1) 



s'^(t)s'M0)=5pas^^ Sit-t'). 



This approximation is valid below the threshold of coil- 
stretched transition, when the polymers do not affect the 
turbulent statistics. This regime will be referred to as the 
"passive" regime. The fourth-rank tensor Spas defined 
by Eq. (3.1) is symmetric with respect to permutations 
of the two first and two last indices, S";-'-' = S";-' ^. 
Incompressibility leads to the restriction Sp^'^*^ = 0. Ho- 
mogeneity implies EfJ^' = EfJ'^. 

At this point we assume also Gaussianity of the turbu- 
lent statistics. Then the tensor Q^ can be found using 
the Furutsu-Novikov decoupling procedure developed in 
[12, 13] for Gaussian processes: 



Qij '^ij^kk^yrkk' 
T ~ "pas -'■^0 



(3.2) 



The cross-correlation tensor Qt is proportional to a 
presently undetermined stress tensor Ho. To find TIq 
one has to substitute Eq. (3.2) into Eq. (2.4) or into its 
formal solution (2.15) and to solve the resulting linear 
system of four equations for the non-zero components of 
Ho: H^^, Uyy, W and H^^ = H?'^. 

The resulting equations are quite cumbersome. The 
basic physical picture simplifies however in two lim- 
iting cases: i) isotropic turbulence. Sec. IIIB, and 
ii) anisotropic turbulence in the limit of strong shear. 

Sec. inc. 



B. Isotropic turbulence 



First we consider the simplest case of isotropic turbu- 
lence. In this case the tensor Hpas has the form 



pas 



pj' -if 
4 V 



, (3.3) 
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where ^ is a constant measuring the level of turbulent 
fluctuations. Using Eq. (3.3) in Eq. (3.2) for one gets 



= s(no5-ino 



(3.4) 



Substituting this relationship into (2.4a) one gets a closed 
equation for Ho: 



\Dt 



1 



So-Uo + Uo-S^ 



+ I SHo + — Heq 

1 S 
2 ■ 



S, 



(3.5a) 



(3.5b) 



Large shear So 3> 1/Tp (Th ^1). In this case the 
threshold velocity gradient is much smaller. In- 
deed: 



I 



4Ife2 



(3.10b) 



Evidently, the strong mean shear decreases the thresh- 
old of the coil-stretch transition very significantly, by a 
factor of Tb^. The important conclusion is that in the 
entire region below the threshold, S < Sc, the renormal- 
ization of the Deborah number can be safely neglected: 



<C 1 



(3.11) 



The stationary solution of Eq. (3.5a) has the form 
Ho = {5 + ?p {So + SI) + 21^So ■ S^) , 



Ho = TV {Ho} . 
In components: 



Ho 



In 

T, 



cq 



(3.6a) 




Ife = T^So<V^ 



(3.6b) 



(3.6c) 



We see that the tensor structure of Ho has the same form 
as in the laminar case, but with an increased relaxation 
frequency given by Eq. (3.5b). 

Taking the trace of Eq. (3.6), one gets the following 
equation for IIq^: 



Txvy 

^^0 



1 



eq ) 



2'De 



(3.7) 

(3.8) 



We see that the effective damping F decreases with in- 
creasing the turbulence level S. At some critical value of 
S = Sc, the effective damping goes to zero, F ^ 0, and 
formally IIq ^ oo. The critical value Sc corresponds to 
the threshold of the coil-stretch transition. Substituting 
Eq. (3.6c) into the threshold condition F = 0, one gets a 
3rd-order algebraic equation for X = ScTp: 



5X^ + 18X^ + 4(3 + AVt^)X = 8 



(3.9) 



This equation has one real root Sc, which we consider in 
two limiting cases: 

• Small shear 5o <C l/rp (Xte ^1). In this case the 
threshold turbulence level Sc is proportional to the 
polymer relaxation frequency 



5td 



1 



(3.10a) 



This means that the structure of the elastic stress ten- 
sor in the passive regime (3.6b) hardly differs from the 
laminar case (2.15). 



The elastic stress tensor in anisotropic 
turbulent field 



Here we consider the case of strong shear, De — > oo, 
having in mind that in the passive regime the tensorial 
structure of Hpas [see Eq. (3.3)] is independent of Tt. In 
this case the leading contributions to each component of 
Ho in Eq. (2.12b) are: 



TTXX 



2rp vc/ Qyy . 



(3.12) 



TTZ2 
^^0 



. 



This means that iovVe^ 1 TVq takes on the form (1.2) 
with 



C: 



■=:vy,xx 



(3.13) 



The threshold condition for H can be found along the 
lines of the preceding subsection. The result is that only 
one component of the tensor Hpas is important: 



:;yy,xx 



1 



2TpIfe2 



(3.14) 



This result coincides to leading order with Eq. (3.10b) in 
which S^lf ^ ^ S according to Eq. (3.3). 

Finally, we point out that the results presented in 
this section remain vaUd for finite, but small, decor- 
relation time Tcor ^ in which S can be defined as 
E =2 f7o2rcor/15, where = ((a"u^)2) = Tr (s • s^). 

Notice that experimental and numerical studies [2, 4, 
16] show that the level of turbulent activity in turbu- 
lent polymeric flows is of the same order as in Newto- 
nian flow at the same conditions. Simple estimations 
for the typical conditions in the MDR regime, when one 
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observes large drag reduction, show that the parameter 
S is far above the threshold value (3.10b). Therefore, 
polymers in the MDR regime cannot be considered as 
passive: there should be some significant correlations be- 
tween polymers and fluid motion that prevent polymers 
from being infinitely extended despite of supercritical 
level of turbulence. The character of these correlations is 
clarified in the following Section. 



IV. ACTIVE REGIME OF THE POLYMER 
STRETCHING 



In this Section we continue the discussion of the case of 
homogeneous turbulence with mean shear, but allow for 
supercritical levels of turbulent fluctuations, at which the 
polymers are strongly stretched. In that case they can 
no longer be considered as a passive field that does not 
affect the turbulent fluctuations. In this active regime 
of the polymer stretching one cannot use Eq. (3.2) for 
the cross-correlation tensor of the turbulent veloc- 
ity field v(r,t) and polymeric stretching Tr{r,t). In the 
present Section we reconsider the v — tt correlations and 
show that in the active regime they are determined by the 
so-called hydroelastic waves which are involved in trans- 
forming turbulent kinetic energy into polymer potential 
energy and vice versa. Therefore in order to find Qt in 
the active regime we need first to study the basic prop- 
erties of the hydroelastic waves themselves. This is done 
in Sects. IV A and IV B. In Sec. IV A we demonstrate 
that the coupled Eqs. (4.3) for v and tt give rise to prop- 
agating hydro-elastic plane waves 



■Ufe , 6fe oc exp[i (fe • r - ujkt) - ^kt] 



with the dispersion law Wfe and damping 7^: 



7fc 



Vik • Ho • fe) 

-(- 
2 lr„ 



(4.1a) 

(4.1b) 
(4.1c) 



For large shear, Xfe ^ 1, these waves (with the exception 
of those propagating exactly in the stream- wise direction 
x) have high quality-factor: uj^^ 7k- The desired cross- 
correlation tensor is defined by the polarization of 
the hydro-elastic waves, that is the subject of Sec. IV B. 
In Sec. IV C we derive the following equation for tensor 

Qt 



— "act ■'■■'■0 



(4.2) 



that is very similar to the corresponding Eq. (3.2) in the 
passive regime. However in Eq. (4.2) the proportionality 
tensor Hact, given by Eq. (4.11a), is very different from 
the corresponding tensor Hpas in the passive regime. 



A. Frequency and damping of hydro-eleistic waves 

The equations (2.9) for the vector v and tensor tt can 
be reformulated in terms of a new vector: 

b = Tp V • TT , (4.3a) 

instead of the tensor tt. The new equations read 

Dv 

— = -So-v + uoAv-Vp + 'ypb + //y, (4.3b) 



Db 1 ^ 

-pr- = b + So-b-Tp^^v +Arb, 

Ut Td 



where 



02 = -V-Ho-V, 

Jsf^ = -V VV, Jsfb = Tp V • 



(4.3c) 



(4.3d) 
(4.3e) 



The linearized version of (4.3) gives rise to hydro-elastic 
waves [9, 14] with an alternating exchange between the 
kinetic energy of the carrier fluid and the potential energy 
of the polymeric subsystem. In the simplest case of space- 
homogeneous turbulent media (no mean shear, 5*0 — 0) 
the homogeneous Eqs. (4.3) have plane wave solutions cf. 
Eq. (4.1). 

When the mean shear exists it appears in Eqs. (4.3) 
with a characteristic frequency 5*0 • In the region of pa- 
rameters that we are interested in, Sq is much smaller 
than the wave frequency Wfc, but much larger than the 
wave damping frequency ^k- This means that the shear 
does not affect the wave character of the motion, but it 
can change the effective damping of the plain wave with 
a given wave vector fe. Indeed, due to the linear inhomo- 
geneity of the mean velocity (constant shear) the wave 
vector becomes time dependent according to 



(4.4) 



This means that in the case ^ jk the shear frequency 
serves as an effective de-correlation frequency instead 
of 7fc: 



Ik ^ 7 = bSo , 



(4.5) 



where b is the dimensionless constant of the order of 
unity. 

We reiterate that the limit Xte » 1 is consistent with 
the frequency of the hydro-elastic waves Wk being much 
larger than their effective damping. 



B. Polarization of hydro-elastic waves 



The formal solution of the linear Eq. (4.3c) for b can 
be written in terms of the Green's function Gpi 

b = -Tp (dp + So ■ dfj V , (4.6a) 



D 1 



(4.6b) 
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Denoting by Fn{t,r) the result of n-fold action on 
the function f{t,r), 



F„(i,r) = G^/(i,r), 



(4.7a) 



one writes 



oo 

Fi {t, r) = j drfit - r, r - So • rr) exp , 



(4.7b) 



By straightforward calculations one can show that Gp 
satisfies the following commutative relationship 

GpV - VGp = • GlV . (4.8) 

Using this relation in Eq. (4.6a) repeatedly we can rewrite 
b as 

(CpV • Ho • V) + So • (g^ V • Ho • V)] v 
(V • Ho • GpV) + (V • So • Ho • G2 v) 
+So • (V • Ho • G2 v) (4.9) 
+2So • (V • So • Ho • G^ v) 



This is an exact relationship for the polarization of the 

hydro-elastic waves in the presence of mean shear, and it 
will be used in the analysis of the cross-correlation tensor 

Qt. 



C. General structure of the cross-correlation tensor 

The cross-correlation tensor Eq. (2.4c) can be 
rewritten in terms of the vector field b as follows: 



Qt = — -{vb + bv) . 



(4.10) 



Substituting the vector b from Eq. (4.9) into Eq. (4.10) 
one finds that the cross-correlation tensor Q^. is pro- 
portional to the elastic stress tensor IIo, according to 
Eq. (4.2), in which the proportionality tensor Hact is ex- 
pressed in terms of the second-order correlation functions 
of the velocity gradients, as follows: 



-■act 



•^ijkl . .^jikl , / -^ijk'l , ■^jik'l\ nk'k 

+ (Si^'^i'''^ + Sf Sf^') (4.11a) 



'^k k 



Here H„ are time-integrated tensors defined as follows: 



(4.11c) 



and r' = r — T So ■ r. 

In order to estimate the relative importance of Si, 
S2, and S3 in Eq. (4.11a) notice that the integrals in 
Eq. (4.11b) arc dominated by contributions from the 
longest hydro-elastic waves in the system. These longest 
fc- vectors have a decorrelation time 7 [which is estimated 
in Eq. (4.5)] , and are characterized by a frequency u that 
we will be specified below. Using these estimates we can 
approximate the time dependence of S'-''^' (r) as follows 



ijkl 



H(0) (4.12a) 



/(r) = cos {u! t) exp (— 7r) 



(4.12b) 



Note that the assumption here is that lo and 7 do not 
depend on the tensor indices. This is a simplifying as- 
sumption that does not carry heavy consequences for the 
qualitative analysis. 

Under these assumptions the tensors H„ in Eq. (4.11a) 
can be estimated as follows 



Z^ijkl J* „ijkl 

fn = 



/i - 



yn-L) 
Re(7 + iujY 
(72 + J^T ' 

/2 ^ -: 



(4.13a) 



(4.13b) 



/3^-^. (4.13c) 



At this point we need to estimate lo. Prom Eq. (4.1b) we 

see that we need to take the largest of component of Ho, 
and smallest available k vector which will be denoted as 
fcmin- Since we expect the xx component of Ho to be the 
largest component we write 

W-Vn^fcrnin. (4.14) 

We should note at this point that homogeneous turbu- 
lence has no inherent minimal k vector, since there is no 

natural scale. In reality there is always an outer scale 
which is determined by external constraints. For fur- 
ther progress in the analysis of the structure of the cross- 
correlation tensor one need to specify the outer scale of 
turbulence. To this end we will consider in the next Sec- 
tion homogeneous turbulence with a constant shear as 
an approximation to wall bounded turbulence of poly- 
meric solution in the region of logarithmic mean velocity 
profile. 



V. POLYMER STRETCHING IN WALL 
BOUNDED TURBULENCE 

In this section we consider turbulent polymeric solu- 
tions in wall bounded flows, and show that the elastic 
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stress tensor takes on the universal form (1.2). This im- 
plies very specific dependence of the components of the 
elastic stress tensor on the distance from the wall. The 
theoretical predictions will be checked against numerical 
simulations and will be shown to be very well corrobo- 
rated. Since we are interested in drag reduction we must 
consider here the active regime when the polymers are 
sufficiently stretched to affect the turbulent field. As 
mentioned before, large drag reduction necessarily im- 
plies Tt ^ 1. For concreteness we restrict ourselves by 
considering the most interesting logarithmic-law region. 
Extension of our results to the entire turbulent boundary 
layer is straightforward. 



A. Cross-correlation tensor in wall bounded 
turbulence 

In this section we consider in more details tensor Qt 
for wall bounded turbulent flows. In this case the outer 
scale of turbulence is estimated as the distance to the 
wall. Therefore in Eq. (4.14) kmin ~ l/y where y is the 
distance to the wall. 

Also we can use the fact that in the turbulent boundary 
layer the mean velocity has a logarithmic profile for both 
Newtonian and viscoclastic flows, (with slopes that differ 
by approximately a factor of five). Therefore the mean 
shear 5o is inversely proportional to the distance to the 
wall y and can be estimated as 

So^—, (5.1) 

y 

where V is the total flux of the mechanical momentum 
(for example, in the channel of half-width L, equal to 
p' L, where p' is the pressure gradient in the streamwise 
direction) . 

Another well established fact is that when the effect of 
drag reduction is large, the momentum flux toward the 
wall is carried mainly by the polymers. Then, Eq. (2.8) 
gives 

nS^~P~(5ot/)^ (5.2) 

To estimate the frequency (4.14) we need to handle 

the other component of the clastic stress tensor. Exam- 
ining Eqs. (2.12b) in the limit Xte ^ 1 we will make the 
assumption that the inequalities (2.13) hold also in the 

strongly active regime. This assiimption will b(^ justified 
self-consistently below. It then follows immediately that 

^o^^s-^ft(5o^/)^ n^^^ng^^^:^. (5.3) 

Now we can estimate the characteristic frequency of 
hydro-elastic waves in Eq. (4.14) as follows: 

w = a V Ve So , (5.4) 



where the dimcnsionlcss parameter a ~ 1. One sees that 
uj indeed is much larger than 7 = 65*0, Eq. (4.5), as we 
expected. 

Next we can continue the analysis of the cross- 
correlation tensor in wall bounded turbulence. Using the 
estimates (4.5), (4.13) and (5.4), we notice that 

Si ~ S0E2 ~ S^VeEs » S^Es . 

Therefore one can neglect terms with S3 on the RHS 
of Eq. (4.11a). Moreover, we can further simplify 
Eq. (4.11a), taking into account that the solution of the 
system of Eqs. (2.4) and (4.11a) preserves the structure 
of polymer stress tensor (2.12b) in which IIq^ ^ ■ 
This allows one to neglect the two last terms in the first 
line of the RHS of Eq. (4.11a), where only S'^y ^ 0. After 
that the cross-correlation tensor Qt, given by Eq. (4.2) 
in terms of tensor Hact) Eq. (4.11a), can be represented 
via the second order tensor 

^ij = fjijMuki (5.5a) 

as follows: 

Q-r = ^2^JHB + B-)-{B-y)Sc-nB-y)] . 

° (5.5b) 
Recall that wc arc looking for the tensorial structure of 
the cross-correlation tensor Qt in order to find the struc- 
ture of IIo. One sees from Eq. (5.5a) that although the 
tensor Ho is strongly anisotropic (its components differ 
by powers of Xte), once contracted with the tensor a*-'''^' 
of a general form, it gives a tensor B with components 
of the same order in Xte. Then Eq. (5.5b) implies that 
the cross-correlation tensor has components that are 
all of the same order in Tk. This result is even stronger 
than the assumption (2.13) in our derivation (that we 
needed to recapture self-consistently). Having done so 
we can conclude that the elastic stress tensor IIo has the 
universal tensorial structure given by Eq. (1.2). 

Armed with this knowledge we observe that the leading 
contribution to B on the RHS of Eq. (5.5a) is given by 

= (j^3,^^ii^^ ^2Ve^ a'^'^^ng^ . (5.6a) 

According to definition (4.11c) the tensor a^i'^^ = 
S'-^'^^(O) can be evaluated as cW^^ /y'^, where y is the 
outer scale of turbulence (and distance to the wall), the 
Reynolds stress tensor W was defined by Eq. (2.7) and 
a new dimensionless constant c is of the order of unity. 
Thus one has 

B'^ = 2 cft2 W'^yy/ ?/2 . (5.6b) 
Now we can write an explicit equation for Q^: 

Q^ = d^Uiy |w-l[(W^.y)x + x(W.y)]| , 

^ (5.7) 

where d = Abc/a^. 
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B. Explicit solution for the polymeric stress tensor 

Having Eq. (5.7) for Qt in terms of XIq and W we 

can find an explicit solution for the mean elastic stress 
tensor Hq in the presence of intensive turbulent velocity 
fluctuations and strong shear. As a first step in Eq. (2.4b) 
for Q we neglect the equilibrium term n^qS/Tp since it is 
expected to be much smaller than Q^, which stems from 
turbulent interactions. Then we substitute Q = Qt into 
Eq. (2.11) [or into Eq. (2.12a)] and solve the resulting 
equations. To leading order in Xfe = SoTp » 1 (in each 
component) this solution takes the form Eq. (1.2) with 



c = w'/wyy . 



(5.8) 



Notice that in our approach the "constant" shear has to 
be understood as a local, y-dependent shear in the tur- 
bulent channel flow, according to Eq. (5.1). Correspond- 
ingly, the Deborah number also becomes ^/-dependent. 



De Ve{y) = So{y)Tp 



\fVi_ 

y 



(5.9) 



Now Eqs. (5.2) and (5.3) provide an explicit dependence 
of the components of tensor Ho on the distance from the 
wall: 











no 




— ) 






y 






— ^h) 




(5.10) 


ttVV 


— >->-0 











At this point we can summarize our procedure as fol- 
lows. First, we assumed that the Q tensor is not strongly 
anisotropic such that weak inequalities (2.13) are valid. 
This allowed us to use the imiversal form (1.2) of ten- 
sor IIo in actual calculations of d the cross-correlation 
function Qt- Then we demonstrated that Qt indeed 
satisfies the required inequalities (2.13). This means that 
Eq. (1.2) presents a self- consistent solution of the exact 
equations in the limit Xte ^ oo. The described procedure 
does not guarantee that Eq. (1.2) is the only solution of 
the problem at hand. However we propose that this solu- 
tion is the realized one, and we will check it next against 
numerical simulations. 



C. Effective polymeric viscosity 

Armed with the structure Eq. (1.2) of the polymeric 
stress tensor, we can rewrite the equation of the mechan- 
ical balance (2.8) as follows: 



v = uoSo-w-y + Tpiiiyso . 



(5.11) 



This means, that the polymeric contribution to the mo- 
mentum flux (last term on the RHS of this equation) can 
be considered as an "effective polymeric viscosity" Vp 



Using Eq. (5.10) one gets universal (rp-independent) lin- 
ear dependence of Vpi 



Up':^VV 



y, 



(5.12b) 



as was suggested in [11] (and sec also [15]). 

The polymeric contribution to the rate of turbulent 
energy dissipation has the form: 



-Tp 



Tr(s.n) = -Tr{Q,}^ 



(5.13) 



see Eq. (A3e) in Appendix A. Using here Eq. (5.7) one 
gets 



b 



(5.14) 



where K = ^Tt{W} is the density of the turbulent ki- 
netic energy. Having in mind that in the MDR regime 
]Y^y ~ K, having the same dependence on the distance 
from the wall, Eq. (5.14) can be rewritten in terms of the 
effective polymeric viscosity fp, given by Eq. (5.12): 



Etp = UpdVul^). 



(5.15) 



fO => t'O + Z^p , t'p = TpUl^ . 



(5.12a) 



where (jVitp) was estimated as K/y"^. 

Notice that the naive estimate for the effective poly- 
meric viscosity is i^p ~ TpHo, that exceeds our re- 
sult (5.12) by a factor of Th"^. The reason for this dif- 
ference is the wave character of the fluid motion; the 
naive result is valid for the estimate of the character- 
istic instantaneous energy flux. As usual in high-quality 
waves or oscillations, the rate of energy exchange between 
subsystems is much larger than the rate of energy dissi- 
pation. 



D. Comparison with DNS data 

The tcnsorial structure of the polymer stress tensor 
was studied in fair detail in DNS for channel flows in 
various papers, see [16-20] and references therein. The 
main problem is that the large 7?e regime at which the 
universal MDR [2] is observed is hardly available. In 
these DNS the maximal available was below 100 at 
the wall, decreasing to about 10 in the turbulent sub- 
layer. For these conditions only up to 50-60 % of the 
total momentum flux is carried by the polymers. Never- 
theless we can compare our analytical results for Ife ^ oo 
with DNS at moderate Tk., at least on a qualitative level. 

The most frequently studied object is Hq = Trjllo}. 
This object is dominated by the stream- wise component 
n§^, see. Fig. 12 in Ref. [16], Fig. 5.3.1 in Ref. [17], Fig. 6 
in Ref. [20], etc. The accepted result is that Hq decreases 
in the turbulent layer. Our Eqs. (5.10) predict Hq oc 1/y 
in this region. This rationalizes the DNS observations in 
Rcfs. [16, 17, 19, 20]. As an example, we present in Fig. 1 
by black squares the DNS data of [19] for the largest. 
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0.0 0.2 0.4 ' 0.6 0.8 1.0 




y 

FIG. 1: Comparison of the DNS data Ref. [19] for the mean 
profiles of the xx and yy components of the elastic stress 
tensor Ilg' with analytical predictions. Black squares: DNS 
data for streamwise diagonal component Hq^, that according 
to our understanding, Eq. (5.10), has to decrease as 1/y with 
the distance to the wall. Black solid line: 1/j/^ dependence. 
Red empty circles: DNS data for the wall-normal component 
Dq", for which we predicted linear increase with y^ in the log- 
law turbulent region. Red dashed line - linear dependence, 
oc 2/+. 



streamwise component Ilg^ and plotted by solid black 
line the expected [see Eq. (5.10)] profile oc l/y"'". In spite 
of the fact that the MDR asymptote with the logarithmic 
mean velocity profile is hardly seen in that DNS, the 
agreement between the DNS data and our prediction is 
obvious. 

The effective polymeric viscosity was measured, for ex- 
ample, in Ref. [16], Fig. 5, and in Ref. [17]. In these Refs. 
fp(j/) was understood as Ap{y)/So{y), where Ap(?/) is 
the Reynolds stress deficit, which is Ilg^ according to 
Eq. (2.8). The observations of Refs. [16, 17] is that 
i^p{y) in the turbulent boundary layer grows linearly with 
the distance from the wall. Recent data of DNS of the 
Navier-Stokes equation with polymeric additives (using 
the FENE-P model) taken from Ref. [19] are presented 
in Fig. 1. The wall- normal component IIq^ is presented 
by the red empty circles. One sees that (in the quite 
narrow) region 20 < < 60 (where the mean velocity 
profile is close of MDR) the profile IIq^ is indeed propor- 
tional to 2/+, as represented by the red dashed line. This 
is in agreement with our result, see Eq. (5.12b) and our 
Ref. [11]. 

Direct DNS of other components of the polymeric 
stress tensor, see, e.g. Fig. 6 in Ref. [19] also agrees with 
the tensorial structure of Ho given by Eq. (5.10). 



VI. SUMMARY AND DISCUSSION 

The aim of this paper was the analysis of the poly- 
mer stress tensor in turbulent flows of dilute polymeric 
solutions. On the face of it this object is very hard to 
pinpoint analytically, being sensitive to complex interac- 
tions between the polymer molecules and the turbulent 
motions. We showed nevertheless that in the limit of 
very high Deborah numbers, Ite 3> 1, this tensor attains 
a universal form. Increasing the complexity of our turbu- 
lent ensemble, from laminar, through homogeneous tur- 
bulence with a constant shear, and ending up with wall 
bounded anisotropic turbulence, we proposed a universal 
form 

Ve{y) \ 
ILo{y) nyy{y) My) 1 . (6.1) 
V C{y)) 

We rewrite this form here to stress that it remains 
unchanged even when the Deborah number, and with it 
the components of the tensor, become space dependent. 
Obviously, this strong result is expected to hold only as 
long as the mean properties, including the mean shear, 
vary in space in a controlled fashion, as for example in 
the logarithmic layer near the wall (be it a Karman or a 
Virk logarithm). 

As an important application of these results we consid- 
ered in Sect. V the important problem of drag reduction 
by polymers in wall bounded flows. A difficult issue that 
caused a substantial confusion is the relation between 
the polymer physics, the effective viscosity that is due to 
polymer stretching, and drag reduction. In recent work 
on drag reduction it was shown that the Virk logarith- 
mic Maximum Drag Reduction asymptote is consistent 
with a linearly increasing (with y) effective viscosity due 
to polymer stretching. This result seemed counter intu- 
itive since numerical simulations indicated that polymer 
stretching is decreasing as a function of y. The present 
results provide a complete understanding of this conun- 
drum. "Polymer stretching" is dominated by Ilg^ since 
it is much larger than all the other components of the 
stress tensor. As shown above, this component is indeed 
decreasing when y increases, cf. Fig. 1. On the other 
hand the effective viscosity is proportional IIq^ , and this 
component is indeed increasing (linearly) with y, cf. Fig. 
1. In fact, drag reduction saturates precisely when Ilg^ 
and Iln^ become of the same order. 
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APPENDIX A: EXACT ENERGY BALANCE 
EQUATIONS 

In this Appendix wc present exact energy balance 
equations that are useful in the analysis of turbulence of 
the polymeric solutions. In the present study we employ 
only one of them, i.e. Eq. (A3e). 

Introduce the mean densities of the turbulent kinetic 
energy and polymeric potential energy Ep 

E^ = ^W, W = Ti{W}, (Ala) 

Ep = ino, no = Tr{no} . (Alb) 

Using Eqs. (2.1) one can derive equations for the bal- 
ance of E-j. and Ep and for their sum: 



DEt 
Dt 



0, 



DE, 



Dt 



— — Sp £p + ^Tp — , 



(A2a) 
(A2b) 



D 



— {E^ + Ep) = 4 +£+-£--£-= . (A2c) 

We denote by e+ and £+ the energy flux from the mean 
flow to the turbulent and polymeric subsystem respec- 
tively: 



4 = -So^y, 

£p = ^>o^i-o , 



(A3a) 
(A3b) 



£^ describes the dissipation of energy in the turbulent 
subsystem, whereas £" is the dissipation in the polymeric 



subsystem due to the relaxation of the stretched polymers 
back to equilibrium: 



= voTr {s ■ s"^) , 

= ^TV(no-neq) 



(A3c) 
(A3d) 



The last term on the RHS of Eqs. (A2a) and (A2b) [that 
is absent in (A2c)] describes the energy exchange between 
the polymeric and turbulent subsystems: 

s-rp = Tr (s • tt) = ^Tr {Q^} . (A3e) 

Using the expression for the momentum flux, we obtain 
an exact balance equation for the total energy 



E = Ey + Erj. + E^ 
that in the stationary state reads: 



p ' 



(A4a) 



SoV = -^Tr (Ho - Heq) + [Sl + Tr (s • s^)) . 

(A4b) 

In Eq. (A4a) E^, is the density of the kinetic energy of the 
mean flow, defined up to an arbitrary constant, depend- 
ing on the choice of the origin of coordinates. The LHS 
of Eq. (A4b) describes the work of external forces needed 
to maintain the constant mean shear. The first term on 
the RHS (oc 1/rp) describes the energy dissipation in the 
polymeric subsystem. The term vi^Sq represents the vis- 
cous dissipation due to the mean shear, while the last 
term on the RHS is responsible for the viscous dissipa- 
tion caused by the turbulent fiuctuations. 
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